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摘要 :西藏 飞 蝗 是 青藏 高 原本 地 特有 物种 ,为 青 委 和 牧草 的 重要 害虫 ,近年 来 危害 范围 有 扩 太 莹 延 趋势 。 研 究 并 明确 西藏 飞 蝗 
在 中 国 的 适 生 区 域 ,对 制定 该 虫 的 早期 监测 预警 及 控制 措施 意义 重大 。 近 年 来 最 大 精 理 论 在 物种 适 生 研 究 领 域 得 到 广泛 应 
用 ,基于 西藏 飞 蝗 的 分 布 信息 和 环境 变量 ,采用 MaxEnt 生态 位 模型 和 AreGIS 对 其 在 趾 国 的 潜在 分 布 区 进行 预测 ,用 ROC 曲线 
检测 模型 精度 、 刀 切 法 (jackknife test) 筛选 主导 环境 变量 及 其 适宜 值 。 两 次 模拟 的 AgC 值 分 别 为 0.996 和 0.993 ,预测 结果 与 实 
际 拟 合 度 很 高 。 西 藏 飞 蝗 在 中 国 的 高 适 生 区 主要 位 于 四 川 的 甘孜 州 ,西藏 的 昌都 地 区 林芝 地 区 山南 地 区 及 拉萨 市 ,中 适 生 区 
则 以 高 适 生 区 为 中 心 向 外 扩散 ,集中 在 青藏 高 原 东 部 地 区 。 海 拔 8 月 份 笠 均 雨 量 :1 月 份 平均 雨量 ,等温 性 12 一 2 月 份 的 平均 
温度 是 影响 西藏 飞 蝗 潜在 分 布 的 主要 环境 变量 。 

关键 词 :西藏 飞 蝗 ;MaxEnt 模型 ; 适 生 分 析 ; 影 响 因子 
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Abstract: Locusta mieratoriadtibetensis Chen, a unique subspecies in the Tibetan Plateau, is a destructive pest of highland 
barley and herbage. The methods of Maximum FEntropy have been deployed for some years to address the problem of species 
abundance distributions. To monitor and control L. migratoria tibetensis, it is necessary to investigate the potential 
distributibfi arealof this pest. In this approach, ecological niche modeling software MaxEnt (the maximum entropy model) ， 
combiiied with_ArcGIS ( Geographic information System), was applied to predict the potential geographic distribution of L. 
niiaratoriatibetensis in China. Bioclimatic dominant factors and the appropriate ranges of their values were also investigated. 
The restilts showed that training data AUC were 0.996 and 0.993 in the two simulations, which indicated a better forecast. 
The highly suitable area for L. migratoria tibetensis was Ganzi in Sichuan Province, and Changdu, Linzhi, Shanan, and 
Lasa of the Tibet Autonomous Region, whereas the moderately suitable area was in the west of Sichuan, east of Xizang, and 


north of Yunnan. The important environmental factors affecting the distribution of L. migratoria tibetensis were altitude, 
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average precipitation in August, average precipitation in January, isothermality, average mean temperature in January ， 


average mean temperature in December, and average mean temperature in February. 
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全 世界 约 有 蝗虫 12000 余 种 , 除 南 极 济 外 各 大 济 均 有 分 布 。 改 凰 (Pocutsa migratoria L.) 是 世界 范围 内 的 
重要 害虫 ,已 知 有 1 个 种 10 个 亚 种 ,其 中 西藏 飞 蝗 ( Locrusta migratoria tibetensis Chen ) 为 青藏 高 原 特 有 贱 分 布 
海拔 最 高 的 亚 种 1 ,主要 分 布 在 西藏 .四 川 及 青海 等 省 区 ,主要 以 青 牺 、 牧 草 及 杂 草 为 食 !*5] 。 据 史料 记载 ， 
西藏 飞 蝗 曾 在 青藏 高 原 多 次 成 灾 ,蝗灾 所 到 之 处 庄稼 颗粒 无 收 。20 世纪 90 年 代 以 来 ,西藏 蝗 在 四 川 甘孜 
州 和 阿坝 州 ,西藏 普兰 .昌都 .阿里 及 青海 玉树 等 地 曾 大 面积 爆发 , 尤 以 2003 一 2006 年 该 虫 发 出 伏 为 严重 , 特 
别 是 雅鲁藏布江 、 雅 敬 江 以 及 金沙 江 等 横断 山脉 河谷 地 带 ,蝗灾 的 爆发 对 当地 农 牧 呈正 常生 产 造 成 极 大 破坏 ， 
导致 难以 估量 的 经 济 损失 !9 。 
ee 预测 物种 适 生 分 布 是 生态 学 的 重要 研究 领域 ,物种 分 布 模型 (species distriutioh models,SDM ) 逐渐 发 展 
并 成 为 研究 物种 适 生性 的 一 个 重要 工具 。 此 类 模型 根据 物种 已 知 的 “存在 2 或 “不 存在 "数据 ,结合 相应 的 
环境 变量 ,利用 一 定 的 数学 理论 计算 物种 在 目标 地 区 的 分 布 概率 , 量 慈 环境 变量 与 物种 分 布 之 间 的 关系 ,寻找 
物种 分 布 的 限制 因子 及 生境 偏好 等 。 目 前 ,物种 分 布 模型 的 常用 领域 主要 为 :预测 物种 潜在 分 布 区 物种 分 布 
二 ” 与 气候 关系 濒危 及 珍贵 物种 适 生 性 预测 和 古生物 地 理 研究 等 方面 “中, 并 起 到 了 良好 的 效果 ''"。 目 前 常用 
-一 ”的 物种 分 布 模型 中 ,最 大 粹 模型 MaxEnt( maximum entropy mg 本 ) 的 模拟 精度 要 高 于 其 他 模型 ,加 之 该 软件 具 
加 有 运行 时 间 短 ,操作 简便 .运行 结果 稳定 和 所 需 样本 量 小 等 特点 在 业内 得 到 了 广泛 的 认可 "中 ,已 被 应 用 于 
多 种 病虫害 适 生 性 的 模拟 预测 ,如 柑橘 黄龙 病菌 ( Candidatus Liberibacter asiaticus ) 、 叶 缘 焦 村 病菌 ( Xylella 
QD fastidiosa) .柑橘 木 虱 ( Diapjporina citri Kuwayama) 和 稳 水 象 甲 (Lissorhoptrus oryzophilus Kuschel) 等 ”| , 目 模拟 
”的 结果 与 物种 实际 分 布 的 吻合 度 高 。 
目前 对 西藏 飞 蝗 的 研究 主要 集中 在 生物 学 特性 ` 环 境 胁迫 .数量 性 状 及 遗传 特征 等 方面 '* ,还 未 发 现 对 
其 生态 位 预测 的 相关 研究 报道 。 本 研究 基于 MaxEnt 模型 ,研究 西藏 飞 凰 分 布 与 环境 变量 之 间 的 联系 ,揭示 影 
响 其 分 布 的 主导 气候 因子 ,分 析 西 藏 展 蝗 在 中 国 的 潜在 分 布 区 , 旨 在 为 制定 合理 的 防 控 措 施 提供 重要 参考 和 
理论 依据 。 


1 材料 与 方法 


1.1 西藏 飞 蝗 的 分 布 数据 与 处 理 

MaxFnt 模型 模拟 需要 物种 分 布 数据 和 环境 数据 ,分 布 数据 即 为 “存在 "点 的 经 纬度 信息 。 西 藏 飞 蝗 的 分 
布 数据 通过 以 下 方式 获得 :(1) 实 地 考察 ,四 川 境内 的 采样 点 多 为 实地 考察 获得 ,调查 点 主要 位 于 西藏 飞 蝗 高 
发 区 的 甘 和 玫 州 和 阿坝 州 ,发现 西藏 飞 蝗 奉 虫 或 成 虫 即 记 为 “存在 ”, GPS 记录 经 纬度 。(2) 查询 物种 分 布 数据 
库 。 本 人 研究 查询 的 数据 库 包 括 “ 国际 农业 与 生物 科学 中 心 (CABI,http ://www.cabi.org/) ”数据 库 “ 全 球 物 种 
多 样 性 信息 库 ( GBIF ,http://www.gbif.org/)”、 中 国 西南 地 区 动物 资源 数据 库 (http://www.swanimal.csdb.cn)、 
教学 标本 资源 共享 平台 (http://mnh.scu.edu.cn/)。(3) 检 索 内 外 公开 发 表 的 西藏 飞 凰 相关 的 期 刊 论文 。 
数据 处 理 : 采 用 组 名 区 分 析 法 对 获得 的 分 布点 进行 校对 、 筛 选 , 排 除 空间 关联 性 较 大 造成 的 过 拟 合 模拟 的 影 
响 ,由 于 环境 变量 空间 分 辨 率 为 2.5 arc-minutes( 约 4.5 km2) , 故 设置 缓冲 半径 为 1.5 km, 当 分 布点 之 间距 离 小 
于 3 km 时 ,只 保留 其 中 一 点 , 共 获 得 分 布点 47 个 ,通过 Google earth 确定 每 个 分 布点 的 经 纬度 。 参 照 MaxEnt 
软件 操作 手册 ,将 经 纬度 坐标 输入 Excel ,保存 格式 为 * .CSV ,字段 包括 物种 名 称 、 经 度 和 纬度 .1 。 由 于 西藏 
飞 蝗 生 活 在 交通 不 便 气候 恶劣 的 高 海拔 地 区 ,所 获得 的 分 布点 很 难 完全 覆盖 其 实际 分 布 。 研 究 表明 ,在 样本 
量 很 小 的 情况 下 ,MaxEnt 模型 较 其 他 模型 预测 效果 更 好 .准确 性 更 高 ,更 适用 于 小 样本 数据 的 预测 。 
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1.2 预测 环境 变量 
1.2.1 环境 变量 数据 来 源 

本 研究 所 用 的 环境 因子 均 来 源 于 Worldclim 气候 数据 库 ,数据 是 以 全 球 的 气象 记录 信息 为 基础 数据 , 整 
合 插值 生成 的 全 球 气候 栅 格 数据 。 利 用 ArcGIS 软件 提取 研究 区 (西藏 四川、 云南 、 贵 州 和 青海 )19 个 生 
物 气候 因子 (具有 较 强 的 生物 学 意义 ,可 反映 温度 与 降水 的 特点 及 季节 性 变化 特征 ) 、 月 平均 降水 量 、 月 平均 
最 高 及 最 低 气温 等 气候 指标 数据 。 另 外 还 选取 了 海拔 数据 作为 地 形 因 子 ( 表 1)。 上 述 数据 的 空间 分 辨 应 5 


arc-minutes( 约 4.5 km”) 。 


表 1 西藏 飞 蝗 潜在 地 理 分 布 的 评价 指标 


Table 1 Climatic variables used for predicting potential geographic distribution of Locusta migratoria tibetenisis 


变量 及 描述 单位 指标 意义 

Variables and description Unit SignificanGe of variaDleg 
Biol = 年 平均 气温 Annual Mean Temperature C 反映 平均 温度 及 其 变异 幅度 
Bio2 = 平均 日 较 差 Mean Diurnal Range( Mean of monthly (max temp-min temp) C 反映 温差 特点 

Bio3 = 等 温 性 Isothermality ( BIO2/BIO7) ( * 100) % 反映 温差 特点 

Bio4 = 温度 季节 性 变化 标准 差 Temperature Seasonality ( standard deviation * 100) @ 反映 平均 温度 及 其 变异 幅度 
Bio5= 最 暖 月 最 高 温度 Max Temperature of Warmest Month C 反映 极端 温度 的 影响 
Bio6= 最 冷 月 最 低温 度 Min Temperature of Coldest Month C 反映 极端 温度 的 影响 
Bio7= 年 均 温 变化 范围 Temperature Annual Range (BIO5-BI06) C 反映 温差 特点 

Bio8= 最 湿 季 度 平 均 温度 Mean Temperature of Wettest Quarter © 反映 水 热 是 否 同步 

Bio9= 最 干 季度 平均 温度 Mean Temperature of Driest Quarter C 反映 水 热 是 否 同步 
Bio10= 最 暧 季度 平均 温度 Mean Temperature of Warmest Quarter C 反映 极端 温度 的 影响 
Bioll= 最 冷 季 度 平均 温度 Mean Temperature of Coldest Quarter EB 反映 极端 温度 的 影响 
Bio12= 年 降水 量 Annual Precipitation mm 反映 降雨 量 及 季节 性 分 布 
Bio13= 最 湿 月 降水 量 Precipitation of Wettest Month mm 反映 极端 水 分 条 件 
Bio14= 最 干 月 降水 量 Precipitation of Driest Month mm 反映 极端 水 分 条 件 
Bio15= 降 水 量变 异 系数 Precipitation Seasonality ( Coefficieiiof Variation) % 反映 降雨 量 及 季节 性 分 布 
Bio16= 最 湿 季度 降水 量 Precipitation of Wettest Quarter mm 反映 极端 水 分 条 件 
Bio17= 最 干 季度 降水 量 Precipitation of DriesteQuarter mm 反映 极端 水 分 条 件 
Bio18= 最 暧 季度 降水 量 Precipitation of W 和 mest Qiter mm 反映 水 热 是 否 同步 
Bio19= 最 冷 季 度 降水 量 PrecipitationvGfticoldest Quarter mm 反映 水 热 是 否 同步 

Alt= 海 拔高 度 Altitude m 反映 海拔 差异 

Tmin= 月 平均 最 低温 度 Minimuin temperattire T 反映 平均 温度 的 影响 
Tmax= 月 平均 最 高 温度 Maximumtemperature C 反映 平均 温度 的 影响 
Tmean= 月 平均 温度 Wean temperature TC 反映 平均 温度 的 影响 
Prec= 月 平均 降雨 最 jPrecipitation mm 反映 平均 雨量 的 影响 


1.2.2 ,影响 西藏 飞 蝗 潜 在 分 布 关键 环境 变量 的 筛选 

影响 物种 分 布 的 环境 因子 多 种 多 样 ,包括 气候 因子 土壤 因子 .植被 分 布 等 ,本 文中 选取 了 worldclim 提供 
的 2 有 个 环境 变量 构建 初始 模型 。 初 始 模 型 构建 过 程 中 选择 MaxEnt 软件 中 的 刀 切 法 (jackknife test ) 检验 来 测 
定 环 境 变量 对 模型 预测 的 贡献 大 小 ,以 此 剔除 对 MaxEnt 模型 预测 结果 贡献 较 小 的 环境 变量 ,对 贡献 较 大 的 关 
键 限 制 因 子 进 行 Spearman 相关 分 析 ,提高 模型 模拟 的 精度 ,在 此 基础 上 重建 西藏 飞 蝗 分 布 的 最 大 焙 模 型 ,并 
对 模拟 结果 进行 准确 性 评价 3] 。 

刀 切 法 检验 使 用 MaxEnt 软件 ,Spearman 相关 分 析 使 用 Spss 分 析 软 件 。 
1.3 MaxFnt 模型 的 操作 与 结果 评价 

MaxEnt 模型 来 源 :选用 基于 最 大 理 论 的 MaxEnt 模型 软件 来 预测 西藏 飞 蝗 在 中 国 的 适 生 区 分 布 ,该 软 
件 可 在 MaxEnt 主页 (http://www. cs. princeton. edu/ ~ schapire/maxent/) 上 免费 获取 "3 ,当前 最 新 版 本 为 
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3.3.3k 版 。 

MaxEnt 模型 的 使 用 :将 上 述 1.1 中 47 个 西藏 飞 蝗 分 布点 数据 和 1.2 中 的 环境 因子 数据 添加 到 软件 中 , 随 
机 选取 75% 的 西藏 飞 蝗 分 布点 作为 训练 集 (training data ) 建立 预测 模型 ,剩余 25% 的 西藏 飞 蝗 分 布点 作为 测 
试 集 (test data) 验证 模型 ,选择 刀 切 法 (jackknife test) 测 定 各 变量 权重 ,选择 创建 环境 变量 响应 曲线 ,其 余 参 数 
均 选择 模型 的 默认 值 '*”|。 

本 人 研究 采用 受 试 者 工作 特征 曲线 ( Receiver Operating Characteristic Curve, ROC 曲线 ) 下 面积 值 , 即 AUC 
值 ( Area Under Curve，AUC ) 为 判 据 来 评价 模型 模拟 结果 。ROC 曲线 ,又 称 为 感受 性 曲线 ,该 曲线 根据 二 系列 
不 同 的 二 分 类 方式 ,分 别 以 假 阳性 率 (1- 特 异 率 ) 和 真 阳 性 率 (1- 遗 漏 率 ) 为 横 纵 坐标 绘制 而 成 ,曲线 让 面积 
值 即 AUC 值 , 取 值 范围 为 [0,1] ,AUC 值 越 接近 1 表明 环境 变量 与 分 布 模型 之 间 相 关 性 越 大 , 锋 测 结果 精度 
越 高 :5 。ROC 曲线 可 使 用 Spss 分 析 软 件 获 得 ,本 研究 中 MaxEnt 模型 可 自动 生成 ROC 曲线 8 并 给 出 相应 的 
AUC 值 。 具 体 评价 标准 站 见 表 2。 


表 2 AUC 评价 标准 表 
Table 2 The evaluation criterion of AUC 


AUC 取 值 范围 评价 标准 AUC 取 值 范围 评价 标准 

Range of AUC values Evaluation criterion Range of AUC value8 Evaluation criterion 
0.5<AUC<0.6 失败 0.8<AUC<09 好 
0.6<AUC<0.7 较 差 0.9<AUE<1.0 极 好 
0.7<AUC<0.8 一 般 


AUC:ROC 曲线 下 面积 值 ， 即 Area Under Curve 


1.4 西藏 飞 蝗 适 生 等 级 分 区 
MaxEnt 模型 输出 结果 为 西藏 飞 蝗 在 世界 范围 内 的 存在 概率 ,数据 为 ASC 本 格式 ,首先 使 用 ArcGIS 的 
ArcToolbox 的 格式 转换 工具 ,将 该 数据 转 为 Rastef 格 式 , 使 该 结果 可 在 ArcGIS 中 显示 ,而 后 使 用 “提取 分 析 ” 
功能 得 到 该 虫 在 中 国 的 存在 概率 分 布 图 :3 WiaxEnt 软件 模拟 输出 的 结果 值 在 0 一 1 之 间 , 值 越 接近 1 表示 
物种 越 可 能 存在 。 参 考 IPCC 报告 1 关于 评估 可 能 性 的 划分 方法 ,结合 西藏 飞 蝗 的 实际 情况 ,利用 
“Reclassify” 功 能 ,划分 分 布 值 等 级 及 相应 分 布 范 围 ,并 使 用 不 同 颜色 表示 ,划分 标准 为 :存在 概率 <0.05 为 不 
适 生 区 ;0.05 三 存在 概率 <0.33 为 低 适 生 区 ;0.33 三 存在 概率 <0.66 为 中 适 生 区 ;存在 概率 三 0.66 为 高 适 生 区 。 


2 研究 结果 


2.1 主导 环境 因子 筛选 及 模型 检验 
2.1.1 影响 西藏 攻 蝗 谐 在 分 布 主导 环境 变量 筛选 

MaxEnt 模型 预测 过 程 中 需 使 用 物种 空间 分 布点 数据 和 环境 变量 数据 ,并 不 是 每 个 环境 变量 对 物种 潜在 
分 布 预 测 都 是 必需 的 ,因此 需要 别 除 对 MaxEnt 模型 预测 结果 贡献 较 小 的 环境 变量 , 簿 选 得 到 贡献 较 大 的 关键 
限制 因 予 8? 。 本 研究 参考 雷军 成 等 “| 方法 ,利用 MaxEnt 软件 本 身 提供 的 模块 , 即 刀 切 法 测定 变量 的 重要 
性 ,通过 比较 每 个 环境 变量 对 模型 结果 的 影响 程度 来 剔除 影响 较 小 的 变量 ,最 终 获 得 对 模型 模拟 贡献 率 较 大 
的 变量 作为 主导 变量 。 本 研究 中 ,依据 MaxEnt 的 运算 结果 ,不 同 环境 变量 对 西藏 飞 凰 潜在 分 布 的 累积 贡献 率 
昂 表 3。 海 拔高 度 为 限制 西藏 飞 凰 分 布 的 主要 变量 , 单 因子 贡献 率 达 到 了 52% ,其 次 ,8 月 份 平均 雨量 、1 月 份 
平均 雨量 ,等温 性 1 月 份 平均 温度 、12 月 份 平均 温度 和 2 月 份 平 均 温 度 对 西藏 飞 蝗 潜在 分 布 的 累积 贡献 率 为 
42.0% 。 

本 研究 使 用 Spearman 相关 系数 法 计算 7 个 环境 因子 之 间 的 相关 性 ,以 消除 共 线性 对 模型 建 模 过 程 和 结 
果 解 释 的 影响 ,环境 因子 之 间 的 相关 系数 绝对 值 大 于 0.8” , 则 认为 两 者 具有 很 强 的 相关 性 。 由 表 4 可 以 
看 出 ,7 个 环境 因子 两 两 之 间 相 关系 数 均 未 超过 0.8 ,因此 本 研究 选 定 上 述 7 种 环境 变量 作为 影响 西藏 飞 蝗 分 
布 的 主导 变量 ,在 此 基础 上 重建 西藏 飞 蝗 在 中 国 的 分 布 的 最 大 炉 模型 ,并 进行 模拟 结果 的 精度 评价 ,以 获得 最 
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优秀 的 预测 模型 。 


表 3 影响 西藏 飞 蝗 分 布 的 环境 变量 贡献 百分率 和 累积 贡献 百分率 


Table 3 The accumulated contribution of each environmental variable to the potential distribution of Locusta migratoria tibetensis defined 


by MaxEnt 
环境 变量 贡献 百分率 累积 贡献 百分率 
Environmental variables Percent contribution/ % Accumulated percent contribution/ % 
海拔 (Altitude，Alt) 5 52 
8 月 份 平均 雨量 ( Average precipitation in August，Prec8 ) 17.3 69.3 
1 月 份 平 均 雨 量 (Average precipitation in January, Precl) 5.9 74.2 
等 温 性 (Isothermality ，Bio3 ) 13 87.2 
1 月 份 平均 温度 ( Average mean temperature in January, Tmean1) 3.4 90.6 
12 月 份 平均 温度 ( Average mean temperature in December, Tmean12) 2.5 93.1 
2 月 份 平均 温度 ( Average mean temperature in February, Tmean2) 1:S 94.6 


累积 贡献 百分率 为 贡献 百分率 之 和 


表 4 环境 因子 之 间 的 Spearman 相关 系数 


Table 4 Pair-wise Spearman’s correlation coefficients of environmental variables 


环境 因子 Alt Bio3 Precl Prec® Tmeanl Tmean2 
Bio3 0.215“ 
Precl 0.134 0.431” 
Prec8 0.15” 0.612” 0.36™ 
Tmeanl 0.017 0.196” 0.115” 0.155 
Tmean2 0.072“ 0.262* 0.142> 0.192“ 0.287 
Tmean12 0.049 0.28“ 0.158 0.196” 0.323 0.023 


Alt :海拔 , Altitude; Bio3 : 等 温 性 ， Isothermality; Precl: 1 月 份 平 均 雨 量 ， Average precipitation in January; Prec8:8 月 份 平均 雨量 ， Average 


precipitation in August; Tmeanl:1 月 份 平 均 温 度 ， Average mean temperatuie in January; Tmean2:2 月 份 平均 温度 , Average mean temperature in 


February;Tmean12 :12 月 份 平 均 温 度 ,Average mean temperéture in Décember; * 表示 在 w= 0.05 水 平 上 显著 相关 


2.1.2 ”ROC 曲线 和 AUC 值 对 模型 准确 性 检验 

图 1 是 MaxEnt 模型 模拟 输出 的 ROC 曲线 图 ,图 中 表明 基于 所 有 环境 变量 和 主导 环境 变量 构建 的 西藏 飞 
凰 在 中 国 的 地 理 分 布 模型 的 AUC 值 分 别 为 0.996 和 0.993 ,依据 表 1 的 AUC 值 评 价 标准 表 , 本 次 构建 模型 的 
预测 准确 性 达到 “ 极 好 的 标准 说 明 此 结果 可 用 于 研究 西藏 飞 蝗 的 适 生 区 划 。 


1.0 上 1.0 
全 0.8 上 0.8 
得 = 
娶 寺 06 上 0.6 
型 
= 
妃子 ”0.4 上 0.4 
时 过 
多 02r 训练 数据 (AUC = 0.996) ”0.2 AUC 平 均值 (AUC = 0.996) 
a 测试 数据 (AUC = 0.995) = 标准 误差 : 士 1 
。 随机 预测 


" 随机 预测 (AUC = 0.5) 


0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 


1- 特 异 率 (小 数 预 测 面积 ) 


1-Specificity (fractional predicted area) 


图 1 模型 适用 性 检验 AUC 值 
Fig.1 AUC value of test on applicability of the MaxEnt model 
A 代表 基于 所 有 环境 变量 计算 的 ROC 曲线 ,B 代表 基于 主导 环境 变量 ROC 曲线 
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2.2 西藏 飞 蝗 在 中 国 的 地 理 分 布 的 预测 

根据 7 个 主导 环境 变量 与 西藏 飞 蝗 的 分 布 数据 构 9 
建 的 最 大 炉 模 型 模拟 结果 ,按照 高 适 生 区 中 适 生 区 低 3NFA < 
适 生 区 和 不 适 生 区 4 个 等 级 标准 划分 ,最 终 获 得 西藏 信 so[”  ， : 
蝗 在 中 国 的 适 生 区 划 图 (图 2) 。 预 测 结果 显示 ,西藏 飞 a 
凰 在 当代 气候 条 件 下 预测 的 适 生 区 范围 为 88.79" 一 ~ a 
113.57?E ,21.52° 一 38.33°*N。 其 中 高 适 生 区 主要 位 于 四 26 上 
川 和 西藏 两 省 区 ,包括 四 川 的 甘孜 州 ,西藏 的 昌都 地 区 、 
林芝 地 区 .山南 地 区 和 拉萨 市 。 中 适 生 区 治 高 适 生 区 有 周 a 
围 分 布 ,主要 集中 青藏 高 原 东 部 地 区 ,包括 四 川西 部 西 适 生 区 划 有 

不 适 生 区 号 低 适 生 区 国 中 适 生 区 el 高 适 生 区 

藏 东部 .云南 北部 等 地 。 

西藏 飞 凰 当代 适 生 区 总 面积 为 170.55x10” km” , 占 2 ”基于 MaxEnt 模型 预测 的 西藏 飞 量 在 中 国 的 适 生 分 析 图 
中 国 国 士 面积 的 17.77% 其 中 ， 高 、 中 低 适 生 区 面 职 Fig. 2 了 Potential suitable sdistributiomzof Locusta migratoria 
分 别 为 22 51x 104 km 人 63.04x 104 km 和 85 01x 104 tibetensis in China based on MaxEnt model 
km2 ,分 别 占 总 适 生 区 面积 的 13.2% 、36.96% 和 49.84%。 四 川 和 西 戏 高 适 生 区 面积 最 天, 分别 为 10.19x10” 
km?* 和 11.05x10* km*( 表 5)。 
2.3 ”西藏 飞 蝗 地 理 分 布 与 环境 变量 之 间 的 关系 
2.3.1 思 切 法 检验 结 日 除 此 变量 a 仅 此 变量 ” 口 所 有 变量 

刀 切 法 检验 结果 可 反映 不 同 环境 变量 对 于 分 布 增 


22” 上 上 


90° 96° 102° 108° 114°E 


三 Alt 
益 的 贡献 大 小 ,该 方法 分 别 计算 “ 仅 此 变量 "“ 除 此 变 pl 
量 " 和 “所 有 变量 "模拟 时 的 训练 得 分 情况 ,判定 标准 凡 Pe 
为 :“ 仅 此 变量 "时 得 分 较 高 ,说 明 该 因子 具有 较 高 的 预 闪 让 Teal? 
测 能 力 , 对 物种 分 布 贡献 较 大 ;“ 除 此 变量 ”模型 训练 。 上 一 一 一 一 一 一 
得 分 能 力 降低 较 多 ,说 明 该 变量 具有 较 多 的 独 等 信息 ， 检疫 但 分 KEiatbejtainbleain 
对 物种 分 布 较为 重要 。 图 3 为 使 用 Wai 软件 自 带 的 
刀 切 法 分 析 各 环境 变量 对 西藏 飞 晶 分布 重 要 性 的 结 ee 


Fig.3 Importance of environmental variables to Locusta 


图 , 申 图 中 可 见 海拔 ( Al) 是 影响 西藏 飞 蝗 分 布 的 关键 
环境 变量 ,其 训练 增益 超过 2. 宇 等温 性 (Bio3 ) 也 是 较为 
重要 的 环境 变量 ,其 单独 训练 得 分 超过 1.5; 另 外 5 种 主导 变量 的 重要 性 排序 为 2 月 份 平均 温度 (Tmean2 ) >1 
月 份 平均 温度 (Timean1) >12 月 份 平 均 温 度 (Tmean12) >1 月 份 平均 雨量 (Precl ) >8 月 份 平均 雨量 (Prec8) ,这 
5 种 环境 变量 的 训练 得 分 也 均 超 过 1.0。 
2.3.2 ”环境 变量 对 MaxEnt 模型 预测 的 影响 

图 4 是 MaxFat 模 型 绘制 的 主导 环境 变量 与 分 布 概率 之 间 的 反馈 曲线 ,该 图 可 反映 不 同 阔 值 下 环境 变量 
的 取 值 范围 。 参 照 1.4 中 的 分 级 方法 ,本文 以 0.33 为 阔 来 划分 适合 西藏 飞 凰 分 布 的 环境 变量 的 范围 。 结 果 显 
示 ., 海拔 的 适宜 范围 为 2870 一 5188 m ,最 适 高 度 为 3941 m ,海拔 2870 一 3941 m 时 ,分 布 概率 随 海拔 升 高 而 增 
加 ,3941 一 5188 m 时 , 随 海 拔 升 高 而 降低 。 等 温 性 的 适宜 值 为 38.5 一 47.2 ,最 适 值 为 42.9 ,在 38.5 一 42.9 , 随 等 
温 性 上 升 分 布 概率 增 大 ,在 42.9 一 47.2 时 , 随 等 温 性 升 高 分 布 概率 降低 。1、2 .12 月 份 平均 温度 和 8 月 份 平 均 
雨量 与 分 布 概率 之 间 的 反馈 曲线 的 变化 趋势 与 上 述 2 种 变量 类 似 ,只 是 在 变化 幅度 及 取 值 范围 存在 差异 , 具 
体 见 表 6。 在 适宜 范围 内 ,各 主导 环境 变量 的 变化 对 西藏 飞 蝗 的 存在 概率 均 有 一 定 的 影响 ,而 取 值 在 适宜 范 
围 之 外 时 对 该 虫 分 布 影响 的 概率 逐渐 减 小 。 


migratoria tibetensis by jackknife analysis 
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表 5 西藏 飞 蝗 在 当代 气候 条 件 下 的 适 生 区 面积 预测 
Table S Predicted areas for Locusta migratoria tibetensis under current climatic condition 
不 庆生 区 低 适 生 区 中 适 生 区 高 适 生 区 
No suitable area Lowly suitable area Moderately suitable area Highly suitable area 
占 该 省 占 该 省 占 该 省 占 该 省 
省 份 i i 积 比例 /% i f 积 比例 /% a i 积 比例 /% a 面积 比例 /% 
站。 预测 面积 面积 比例 /% 预测 面积 。 西 积 比例 /% 预测 面积 。 。 宣 积 比例 /% 预测 面积 /面积 比例 /% 
Province 人 Account for 站 这 Account for Account for 有 Account for 
X10”km- . x10” km . x10” km . x10” km’” . 
1 the proportion of the proportion of the proportion of the proportioiy of 
Predicted area Predicted area Predicted area Predicted area 
the area of the the area of the the area of the the area_of the 
province province province povince 
川 6.15 13.51 11.03 24.23 18.14 39.86 10.19 22.40 
西藏 30.92 31.45 25.12 25:55 31.22 31.76 11.05 11.24 
云南 533 15.57 19.41 56.73 8.22 24.03 1.26 3.67 
贵州 9.78 61.35 5.50 34.51 0.66 4.14 0.00 0.00 
青海 42.58 59.70 23.94 33.57 4.80 6.73 0:01 0.01 
1.0 Alt 
0.5 
0 
-74 6022 
1.0 Bio3 1.0 Precl 1.0 Prec8 
(3 0.5 0.5 
0 0 0 
9 93 0 756 0 1095 
Tmeanl Tmean12 Tmean2 
1.0 1.0 1.0 
Us 0.5 5 
0 0 0 
-508 330 -490 329 -470 323 


图 4 MaxEnt 模型 中 研究 物种 对 各 环境 变量 的 反馈 曲线 


Fig.4 Response of curves of environmental variables in MaxEnt models 


表 6 西藏 飞 蝗 潜在 分 布 对 应 的 环境 变量 适宜 范围 表 
Table 6/ The suitable range of dominant environmental variables affecting the potential distribution of Locusta migratoria tibetensis 

环境 变量 适宜 范围 最 适 值 
Environmental Variables Suitable range Optimum value 
海拔 Altitudé 2870 一 5188 3941 

等 温 性 Bothermality 38.5 一 47.2 42.9 

1 月 份 平均 雨量 Average precipitation in January 0 一 9.35 3.56 

8 月 份 平均 雨量 Average precipitation in August 79.7 一 296 98.2 

1 月 份 平 均 温 度 Average mean temperature in January -13.3 一 6.73 -3.1 

12 月 份 平 均 温 度 Average mean temperature in December -12.12 一 7.13 -2.13 

2 月 份 平均 温度 Average mean temperature in February —11.18—8.78 -1.51 
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研究 表明 1 ,环境 变量 的 选择 会 对 生态 位 模型 的 预测 结果 造成 一 定 的 影响 。 应 用 MaxEnt 模型 预测 物种 
地 理 分 布 的 很 多 研究 ,采用 世界 气候 数据 库 (http://www.worldclim.org) 的 Bioclim 数据 集中 的 19 个 生物 气候 
因子 作为 全 部 或 主要 环境 变量 。19 个 环境 变量 是 以 温度 和 雨量 为 基础 数据 ,根据 不 同 需求 计算 演 生 而 来 ， 
此 这 些 变量 之 间 存 在 不 可 避免 的 自 相 关 及 多 重 线性 重复 等 问题 :“”; 。 研 究 业已 证 实 ,这 些 高 度 相关 的 变量 
在 模型 预测 过 程 中 会 引入 宛 余 信息 ,影响 预测 结果 。 在 建 模 过 程 中 ,为 避免 上 述 问题 ,应 首先 对 环境 变量 进行 
相关 分 析 和 有 效 筛 选 。 本 研究 中 ,参考 孙 敬 松 等 ' 苦 史 方法 ,以 各 因子 对 物种 分 布 的 贡献 率 来 评价 变量 的 重要 
性 ,有 效 剔 除了 贡献 率 较 小 的 环境 变量 并 对 第 选 的 变量 进行 了 相关 性 分 析 , 获得 了 7 种 主导 环境 谈 量 并 对 模 
型 进行 重建 ,减少 了 宛 余 信息 对 模拟 结果 的 影响 ,提高 了 预测 结果 的 准确 性 。 

目前 对 模型 精度 评价 应 用 最 广泛 的 方法 为 ROC 曲线 法 ( 即 AUC 法 )5 ,由 笠 AUC 不 巡 诊 断 闻 值 的 影 
响 ,并 可 提供 所 有 阅 值 范围 上 的 性 能 评价 结果 ,因此 目前 被 公认 为 生态 位 模型 评价 领域 的 最 佳 评价 指标 。 
AUC 的 取 值 范围 在 0.5 一 1, 值 越 接近 于 1 模型 精度 越 高 。 本 研究 中 基于 主导 环境 变量 的 模拟 的 训练 集 AUC 
平均 值 为 0.993( 非常 接近 1) ,预测 结果 达到 “ 极 好 ”水 平 ,说 明 此 次 模型 预测 的 地 理 分 布 与 西藏 飞 蝗 实际 分 布 
拟 合 度 较 高 。 另 外 ,本 研究 利用 地 理 信息 系统 软件 ArcGIS 对 MaxEh 输 出 的 栅 格 文 翌 进行 后 期 处 理 , 使 目标 
物种 的 分 布 数据 与 环境 变量 数据 在 栅 格 单元 上 相对 应 ,有效 的 降低 了 系统 误差 ,进一步 提高 了 数据 的 准确 性 。 

本 研究 预测 结果 显示 ,西藏 飞 蝗 在 中 国 的 适 生 范围 为 88.79° 一 壬 3.57°E,21.52° 一 38.33°N, 总 面积 达 
170.55x10” km 。 其 中 高 适 生 区 主要 位 于 四 川 和 西藏 两 省 区 46 审 适 引 区 则 主要 集中 青藏 高 原 东 部 地 区 ,包括 
四 川西 部 .西藏 东部 云南 北部 等 地 。 上 述 区 域 地 理 上 包 伟 对 四 川 的 多 数 草 原 (若尔盖 草原 、 红 原 草原 、 塔 公 
草原 ,龙灯 草原 和 毛 干草 原 等 ) 以 及 西藏 雅鲁藏布江 及 雅 芍 江河 谷 。 李 庆 等 5 研究 发 现 , 西 藏 飞 蝗 在 草地 环 
境 下 产 卵 块 数 最 多 .草地 环境 更 适合 其 生活 ,西藏 改 蝗 可 对 草地 造成 严重 和 危害。 因此 西藏 飞 蝗 在 上 述 地 区 的 
危险 系数 极 高 ,有 关 部 门 应 引起 高 度 重 视 , 严 防 西藏 飞 峰 聚集 成 灾 。 封 传 红 等 "5 建立 了 有 效 积温 与 地 理 位 置 
的 模型 ,获得 了 西藏 飞 蝗 的 潜在 分 布 面积 ;结果 表明 ,1998 年 西藏 飞 粕 分 布 面积 最 大 ,为 14.3x104 km2。 本 文 
中 高 适 生 区 面积 为 22.51x10” km ,中 适 生 区 面积 则 达 63.04x10” km ,与 封 传 红 的 研究 结果 相差 较 大 ,这 可 能 
是 因为 使 用 的 模型 及 变量 不 同 所 导致 。 

物种 -环境 关系 是 研究 物种 生态 学 需求 和 空间 分 布 的 重要 方面 ,本 研究 分 析 了 西藏 飞 蝗 存在 概率 与 主导 
环境 变量 之 间 的 关系 ,获得 光 相关 的 反馈 曲线 。 结 果 表明 ,西藏 飞 蝗 的 存在 概率 随 着 主导 环境 变量 (海拔 .等 
温 性 1 月 份 及 8 月 份 平均 雨量 .月份 2 月 份 及 12 月 份 平均 温度 ) 的 变化 而 改变 。7 个 主导 变量 中 ,海拔 和 
等 温 性 是 影响 该 虫 分 布 最 为 重要 的 环境 变量 ,表明 西藏 飞 蝗 扩 散 和 繁殖 过 程 中 ,一 方面 受 当地 海拔 范围 的 制 
约 , 另 一 方面 也 受 当 地 等 温 性 的 影响 。 本 研究 中 西藏 飞 蝗 适宜 的 海拔 范围 为 2870 一 5188m , 最 适 为 3941 m。 
目前 已 知 西藏 飞 帆 的 分 布 范围 仅 限 于 青藏 高 原 地 区 , 王 思 忠 等 '5 虽 曾 在 四 川 雅安 (海拔 560 m 左右 ) 室内 成 
功 饲养 繁殖 ,但 死亡 率 高 , 且 众 多 研究 证 明 至 今 未 发 现 该 虫 在 低 海拔 地 区 分 布 ,这 与 本 实验 的 结果 相 吻 合 。 等 
温 性 网 映 了 温度 变 兹 的 迟早 及 幅度 ,如 当年 气温 转变 迟缓 , 则 利于 西藏 飞 凰 的 豚 胎 发 育 ,西藏 飞 蝗 存 在 概率 较 
高 的 等 温 性 指标 在 42.9 左右 ,这 表明 该 虫 在 温差 较 大 情况 下 存在 概率 较 高 ,这 与 上 述 的 研究 结果 相近 。 李 庆 
等 研究 了 西藏 飞 蝗 各 发 育 阶 段 的 耐寒 性 ,结果 表明 卵 期 的 过 冷却 点 为 -22.02% ,4 龄 蝗 晴 的 过 冷却 点 为 
-6.46% ,西藏 飞 蝗 以 卵 越冬 。 西 藏 飞 蝗 蝗 卵 越冬 期 长 达 200 d, 冬 季 低 温 会 影响 其 胚胎 发 育 及 孵化 率 。 封 传 
红 等 ' 生 研究 发 现 ,温度 为 限制 西藏 飞 蝗 分 布 的 主要 限制 因子 ,夏季 高 温 不 是 限制 因素 ,冬季 低温 可 影响 西藏 
飞 凰 的 安全 越冬 ,这 些 都 说 明 冬 季 低 温 对 蝗 卵 的 越冬 至 关 重要 ,本 研究 认为 西藏 飞 凰 的 存在 概率 与 越冬 期 12 
月 及 次 年 1.2 月 平均 温度 有 关 , 在 较 低 的 温度 条 件 下 ,该 虫 仍 有 较 高 的 存在 概率 ,但 过 低 的 温度 将 导致 其 存在 
概率 下 降 ;8 月 份 平 均 雨 量 也 是 影响 西藏 飞 蝗 分 布 的 重要 因素 ,最 适 值 为 98.2 mm。 研 究 发 现 ,西藏 飞 蝗 与 飞 
蝗 其 他 亚 种 类 似 ,主要 林 居 在 较 低 湿 的 禾 本 科 杂 草地 带 , 该 虫 一 般 在 8 月 开始 交配 、 产 卵 ,适宜 的 湿度 有 利于 
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该 虫 顺利 的 完成 这 些 生命 活动 。 此 处 需要 注意 ,物种 存在 概率 对 环境 变量 的 反馈 曲线 反映 的 是 单一 环境 变量 
的 作用 ,但 昆虫 的 生命 活动 受到 各 种 环境 变量 (包括 气候 因子 、 寄 主 条 件 . 天 敌 种 类 植被 条 件 等 ) 的 综合 影 
响 , 因 此 此 结果 可 作为 判断 西藏 飞 蝗 与 环境 变量 之 间 关系 的 参考 ,但 不 能 完全 准确 解释 二 者 之 间 的 关系 。 

尽管 MaxEnt 模型 预测 具有 操作 简便 样本 需求 量 小 和 预测 精度 高 等 优势 ,但 同 其 他 生态 位 预测 模型 类 似 
不 可 避免 的 存在 一 定 的 局 限 性 。(1) 用 MaxEnt 模型 进行 预测 时 所 涉及 的 环境 变量 除 海拔 外 ,其 余 变 量 均 为 
气候 变量 ,19 个 生物 气候 变量 为 西藏 飞 凰 实际 分 布地 的 气候 极 值 , 即 最 大 值 和 最 小 值 ,因而 ,MaxEnt 模型 所 显 
示 的 是 物种 分 布 的 最 大 可 能 性 ,无 法 准备 表达 物种 在 现实 分 布 的 主要 地 区 。 上 述 预测 结果 更 侧重 于 理解 和 展 
现 西藏 飞 凰 潜在 地 理 分 布 以 及 揭示 适 于 该 物种 分 布 的 气候 特征 。(2) 基础 生态 位 是 一 种 理想 生态 位 ,日 其 六 
物种 在 理想 的 生存 条 件 而 不 存在 其 他 物种 竞争 时 所 占据 的 最 大 生态 位 ,该 理论 只 需 考 虑 非 生物 因素 的 影 
响 '”| 。MaxEnt 模型 对 适 生 区 的 预测 就 是 基于 物种 对 基础 生态 位 的 需求 ,但 在 实际 的 生存 环境 中 性 物 因 子 
(如 物种 之 间 的 相互 作用 、 植 被 类 型 .地貌 特征 、 物 种 自身 扩散 能 力 和 土壤 类 型 等 疗 样 会 对 预测 物种 的 潜在 
分 布 产生 重要 影响 '*。 例 如 , 李 庆 等 站 研究 了 植被 及 土壤 对 西藏 飞 蝗 产 卵 及 分 布 的 影响 凡 果 表明 ,植被 类 
型 .土壤 含 盐 量 .含水 量 和 PH 值 能 显著 影响 西藏 飞 蝗 的 产 卵 及 分 布 。 基 于 上 述 原 因 , 可 以 推测 , MaxEnt 模型 
预测 的 生态 位 比 西藏 飞 蝗 所 占据 的 实际 生态 位 要 宽 。 对 此 ,在 下 一 步 工 作 中 , 除 券 虑 气候 因素 影响 外 ,还 应 注 
重 考量 物种 间 相 互 作用 等 生物 因子 的 可 信 表 达 , 以 改善 模型 的 预测 效果 。 | (3 ) 过 去 39 秆 来 , 随 着 温室 气体 的 
不 断 排放 和 保持 ,全 球 气 候 正 异常 变 暧 ,而 气候 变化 会 引起 物种 生长 及 分 布 梅 局 的 变化 '”” 。 本 研究 中 所 使 
用 的 环境 变量 来 源 于 世界 气候 数据 库 worldclim ,时 间 跨 度 为 195$0 一 2000s 缺失 了 近期 10 多 年 的 气候 数据 。 
因此 ,在 以 后 的 工作 中 ,应 充实 缺失 的 数据 ,使 预测 结果 更 为 准备 可 靠 。 
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